Exploring Topological Phases With Quantum Walks 
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The quantum walk was originally proposed as a quantum mechanical analogue of the classical ran- 
dom walk, and has since become a powerful tool in quantum information science. In this paper, we 
show that discrete time quantum walks provide a versatile platform for studying topological phases, 
which are currently the subject of intense theoretical and experimental investigation. In particular, 
we demonstrate that recent experimental realizations of quantum walks simulate a non-trivial one 
dimensional topological phase. With simple modifications, the quantum walk can be engineered 
to realize all of the topological phases which have been classified in one and two dimensions. We 
further discuss the existence of robust edge modes at phase boundaries, which provide experimental 
signatures for the non-trivial topological character of the system. 
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Quantum walks, the quantum analogues of classical 
random walksW, form the basis of efficient quantum 
algorithms 0, |J and provide a universal platform for 
quantum computation [4]. Much like their classical coun- 
terparts, quantum walks can be used to model a wide va- 
riety of physical processes including photosynthesis [^Q, 
quantum diffusion optical/spin pumping and vor- 
tex transport 0, and electrical breakdown[9i, J^. Mo- 
tivated by the prospect of such an array of applica- 
tions, several groups have recently realized quantum 
walks in experiments using ultracold atoms in optical 
lattices 11 1, trapped ions[12|. photons 13 1, and nuclear 
magnetic resonancefl^. These systems offer the possi- 
bility to study quantum dynamics of single or many par- 
ticles in a precisely controlled experimental setting. 

Here we show that quantum walks can be used to ex- 
plore dynamics in a wide range of topological phases [isl- 
17|. Interest in topological phases was first sparke d by 
the discovery of the integer quantized Hall effect [TtI. llSlj. 
and has r apidly increased in recent years following the 
prediction 1^ 21 1 and experimental realization (22i [isj of 
a new class of materials called "topological insulators." 
Unlike more familiar states of matter such as the fer- 
romagnetic and superconducting phases, which break 
SU(2) (spin-rotation) and U(l) (gauge) symmetries, re- 
spectively, topological phases do not break any symme- 
tries and cannot be described by any local order param- 
eters. Rather, these phases are described by topolog- 
ical invariants which characterize the global structures 
of their ground state wavefunctions. Topological phases 
are known to host a variety of exotic phenomena such as 
fractional charges and magnetic monopoles 24, 25| . 

The class of topological phases which can be realized in 
a system of non-interacting particles is determined by the 
dimensionality of the system and the underlying symme- 
tries of its Hamiltonian. Figure [1] shows the ten classes 
of topological phases which can arise in one dimensional 
(ID) and two dimensional (2D) systems with and without 
time-reversal symmetry (TRS) and particle-hole symme- 
try (PHS) (see Refs.0, [53] and discussion below). If 
both symmetries are absent in ID, the possibility of a 
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FIG. 1: Classification of topological phases by symmetry for 
one (ID) and two (2D) dimensions, adapted from Ref . [2^. 
[27( 1. Discrete time quantum walks can naturally realize all ten 
classes of nontrivial topological phases in ID and 2D. Time 
reversal symmetry (TRS) and particle-hole symmetry (PHS) 
are defined by the existence of antiunitary operators T and V 
satisfying Eqs.([7| and ([8]), and may be absent, or present with 
7^ = ±1 {V^ = ±1). In the absence of both TRS and PHS, a 
distinct "chiral" symmetry with a unitary F satisfying Eq.(l9]) 
may be found. In each case, the symmetry-allowed phases 
are classified by an integer (Z) or binary {Z2) topological 
invariant. Classes containing the Su-Schrieffer-Heeger model 
(SSH)[l|], integer quantized HaU (IQH) [13, [H, and quantum 
spin-Hall (QSH)[13-[2^ phases are indicated. 



distinct "chiral" symmetry creates an additional class of 
topological phases. Within each class, the allowed phases 
are characterized by either an integer (Z) or binary (Z2) 
topological invariant. 

In this paper we investigate the topological phases of 
Fig H] in discrete-time quantum walks (DTQWs). In a 
DTQW, a walker with a two-fold internal "spin" degree of 
freedom is made to hop between adjacent sites of a lattice 
through a series of unitary operations. We discuss how 
the DTQW protocol can be engineered to selectively sat- 
isfy time-reversal and particle-hole symmetries, and show 
that DTQWs can realize all of the classes of topological 
phases in ID and 2D. In particular, we show that the 
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FIG. 2: (a) One dimensional discrete time quantum walk 
(DTQW) protocol. First, the walker's internal "spin" is ro- 
tated through an angle d about the t/-axis. Then, the walker 
is coherently translated by one lattice site to the right (left) if 
its spin is up (down). The quantum walk is produced by re- 
peatedly applying this combined "step" operation, (b) Band 
structure of ID DTQW with d — -k 12. The spinor eigenstates 
at each momentum k are directed along the unit vector ng (fc) 
[see Eq.((6|], as represented on the Bloch sphere in panel (c). 
Corresponding points in panels (b) and (c) are indicated by 
the colored markers. For any Q 7^ 0, 27r, ne(fc) winds around 
the origin once as k traverses the Brillouin zone (black arrow). 
The vector Ae is normal to the plane containing ne(fc), and 
defines the axis for the chiral symmetry Fg, see Eq.([9ll. 



DTQWs demonstrated in recent experiments have 
ah'eady realized a non-trivial one dimensional topological 
phase, which is analogous to that of the Su-SchriefFer- 
Heeger (SSH) model of polyacetylene[15j (see FigH]). 

The non-trivial topological properties of the systems 
classified in FigH] are manifested in the presence of ro- 
bust edge states at phase boundaries, i.e. zero energy 
bound states[l^ and gapless edge modesji^ in ID and 
2D systems, respectively. We propose a scheme to iden- 
tify the presence of topological phases through the ob- 
servation of edge modes at an interface between regions 
where different DTQW protocols are applied. 



TOPOLOGICAL PHASES IN ID 

The one dimensional DTQW protocol employed in 
recent experiments Ill4l4| is depicted schematically in 
FigHh- The basis states of the system are described in 
terms of the position of the "walker," defined on integer 
lattice sites x, and its internal "spin" state which can 
be either up (f) or down (J,). The quantum evolution is 
produced by repeatedly applying a unitary operation 



u{e) = TR{e) 



(1) 



that defines one "step" of the quantum walk. Each step 
consists of a spin rotation i?(0), followed by a coherent 



spin-dependent translation 
T = Y,[\^ + I tXt I + |x - ® I ;)(; I] (2) 

X 

that shifts the walker to the right (left) by one lattice 
site if its spin is up (down). This step protocol is a uni- 
tary generalization of the classical process in which a ran- 
dom walker hops left or right according to the outcome 
of a stochastic "coin-fiip." Here, as in the experiments 
of Refs.jTl. 12 1, we consider the case where R{0) corre- 
sponds to a spin rotation around the y-axis through an 
angle 9, 



m 



cos (61/2) -sin (61/2) \ 
sin (61/2) cos(e'/2) ) 



(3) 



Although the step protocol is defined explicitly in 
terms of the discrete unitary operations T and R{9), the 
net evolution over one step is equivalent to that gener- 
ated by a time-independent effective Hamiltonian H{9) 
over the step-time 5t, 



U{9) 



'iH(e)st 



1. 



(4) 



The evolution operator for N steps is given by (9) = 
^-tH(9)NSt_ rpj^^g^ DTQW provides a stroboscopic 
simulation of the evolution generated by H{9) at the dis- 
crete times N5t. Below we take units in which 5t = \. 

The DTQW protocol described above is translationally 
invariant. The evolution operator U (9) and the Hamilto- 
nian H{9) are thus diagonalized down to 2 x 2 blocks in 
the basis of Fourier modes, \k)(^\a) = -^=J2x^~^''^\^)^ 



|(t), with -TT < k < TT. For the choice of R{9) in Eq.(|31), 
H{9) can be written as 



Hi9) 



dk[Eg{k)ne{k) ■ (T](g)\k){k\, (5) 



where <t = (a^ , cr.y ,172) is the vector of Pauli matrices and 
the unit vector ng{k) = {nx,ny,nz) defines the quanti- 
zation axis for the spinor eigenstates at each momen- 
tum k. Because the evolution is prescribed stroboscop- 
ically at unit intervals, the eigenvalues ±Ee{k) of H{9) 
are only determined up to integer multiples of 27r. The 
corresponding band structure is thus a "quasi-energy" 
spectrum, with 27r periodicity in energy. For 9 ^ or 
27r, explicit expressions for Eg{k) and ng(k) are given by 
cos Eg{k) — cos(f?/2) cosfc and 



(sin((?/2) sinfc, sin(6'/2) cos fc, — cos(6'/2) sin fc) 
sm Eg{k) 

(6) 

A typical band structure ztEg{k) is shown in Figl^Jj. 
Note that for 0* = or 27r, the spectrum of H{9^,) is 
gapless and n^^ (fc,) is ill-defined for fc, = 0, tt. 

Hamiltonians of the form ([5]) can support topological 
phases if they possess certain symmetries, as indicated 
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in Fig[T] The time-reversal and particle-hole symmetries 
of this table are defined by the existence of antiunitary 
operators T and V satisfying 



THT-^ 



-H. 



(7) 
(8) 



The Hamiltonian H{9) given by Eqs.® and (O possesses 
PHS (HI) with V = K, where K is the complex conjuga- 
tion operator. To see this, note that the evolution opera- 
tor U (9) given by Eqs.([l][3l) is real, and thus invariant un- 
der K. Along with Eq.gl), this implies H*{0) = -H{9), 
which satisfies Eq. ([S]) with V = K. In addition, using 
Eq. ([5]), it is straightforward to check that H{d) possesses 
a unitary "chiral" symmetry of the form 



Tg'H{e)Te = -H{6), 



(9) 



with Te = g-i^Ae,-<T/2^ where = 
(cos(6'/2), 0, sin(0/2)) is perpendicular to ng[k) for 
all k. The presence of both PHS ([8]) and chiral symme- 
try ([9]) guarantees that H{0) is invariant under TRS ([7]) 
with r = TgV, see Refs.jll |27|. 

The symmetry classes identified in FiglT] are distin- 
guished by whether the relevant symmetry operators T 
and V square to 1 or —1. Because here both = 1 
and V'^ = 1, H{0) belongs to the class of Hamiltonians 
labeled "SSH." The corresponding integer-valued topo- 
logical invariant Z has a simple geometrical interpreta- 
tion. Chiral symmetry © constrains n0(fc) to lie on a 
plane which is perpendicular to Ag, and which contains 
the origin (see Fig. [21:). Thus, H{6) can be characterized 
by the number of times ne(fc) winds around the origin as 
k runs from — tt to tt. Since the winding number of Ta.e{k) 
given by Eq.® is 1 for all^^^O, 27r, the DTQWs imple- 
mented in experiments 
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simulate the Z = 1 SSH 

topological phase. 

The non-trivial topological character of the system can 
be revealed at a boundary between topologically distinct 
phases. To open the possibility to create such a bound- 
ary, we introduce the "split-step" DTQW protocol shown 
in FigHi. Starting from the DTQW defined by Eq.©, 
we split the translations of the spin-up and spin-down 
components, and insert an additional spin-rotation R{d2) 
around the j/-axis in between: 



(10) 



where 1^(4-) shifts the walker to the right (left) by one 
lattice site if its spin is up (down) . 

The split-step protocol defines a family of effective 
Hamiltonians Hss(Qi,Q2) parametrized by the two spin- 
rotation angles Oi and O2 ■ This family realizes both Z = 
and Z = 1 SSH topological phases as displayed in FiglSb, 
with chiral symmetry (|9]) given by Te^fi^ = V — K, 
and T = Tg-^V (see Methods). Gapped phases with wind- 
ing numbers Z = and Z — 1 are separated by phase 



transition lines where the quasi-energy gap closes at ei- 
ther E = or E = ±77, as indicated in the figure. 

We propose to create a phase boundary in the DTQW 
by replacing the second (spatially-uniform) spin rota- 
tion R{92) of Eq. pU)) with a site-dependent spin rotation 
R[02{x)], which rotates the walker's spin through an an- 
gle 92{x) about the y-axis at each site x. Specifically, we 
consider the situation where 6*2 (x) — ?> 62- for x ^ and 
changes monotonically to 6*2 (x) — > 62+ for x ^ (see 
FiglSt). Although this protocol is not translationally- 
invariant, symmetries ([TUn]) are preserved. In particular, 
the system retains the chiral symmetry under Tg^ for ar- 
bitrary 02 (x) as long as 0i remains uniform. 

When the rotation angles (01,6*2+) and (01,02-) are 
chosen to realize topologically distinct phases with Z = 
and Z = I'm the regions x ^ and x ^ 0, a bound state 
with energy or tt exists near the phase boundary x = 
[l6j . The existence of such a bound state is guaranteed 
by topology, and does not depend on the details of the 
boundary. The bound state can be probed by initializing 
the walker at x = as demonstrated in FiglSjl. Be- 
cause this initial state has a non-vanishing overlap with 
the bound state, part of the walker's wavepacket will re- 
main localized near x = 0. On the other hand, if the 
pairs (01, 02+), (01, 02-) are chosen to lie within the same 
"diamond" -shaped region of FiglSb, then the system can 
be made spatially uniform through a continuous defor- 
mation of the Hamiltonian without closing either gap at 
E = Q 01 E = -K. In this case, there are no topolog- 
ically protected modes localized at the boundary. For 
monotonic 02 (x), this guarantees that the system does 
not support any bound states, and the probability to 
find the walker at x = decays to zero with an increas- 
ing number of DTQW steps (see FiglSjs). 

With further modifications to the DTQW protocol, 
each of the topological classes in ID given in FiglUcan be 
realized (see Supplementary Material). In addition, as we 
will now discuss, a straightforward extension of the pro- 
tocol to a higher-dimensional lattice allows the DTQW 
to simulate topological phases in two dimensions. 



TOPOLOGICAL PHASES IN 2D 

To begin, we consider a family of 2D quantum walks in 
which the walker possesses two internal states as in the 
ID DTQWs above. Non-trivial topological phases can 
be realized in a variety of 2D lattice geometries. Here 
we consider the case of a triangular lattice, and discuss 
equivalent square lattice realizations in the Supplemen- 
tary Material. One step of the quantum walk is defined 
by the unitary operation 

(72d(01, 02) = r3i?(0i)r2i?(02)Tii?(0i), (11) 

where Ti {i — 1,2,3) translates the walker with spin up 
(down) in the +(— )vi direction, with {v^} defined in 
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FIG. 3: (a) One-dimensional "split-step" DTQW protocol, see Ea. (|10p . (b) Winding number associated with the split-step 
DTQW as a function of the spin-rotation angles 9i and 62. Topologically distinct gapped phases are separated by phase 
transition lines where a gap closes at either E = ot E = n. (c) Phase boundary in the spatially inhomogeneous split-step 
DTQW. In the second rotation stage of Ea.(|10[l. the walker's spin is rotated by an angle 02 (a;) = |(02- + 6*2+) -I- — 
02-) tanh(x/3). (d), (e) Dynamics of the spatially inhomogeneous split-step walk, with the walker initialized with spin up at 
X = 0. In both panels, we take 81 — —77/2 and 82- — Sn/A, corresponding to winding number in the region x <^ (white dot 
in panel b). In (d) we create a phase boundary by taking 82+ = 7r/4, which gives winding number 1 for j: S> (see blue dot 
in panel b). After many steps, the probability to find the walker near a; = remains large, indicating the existence of at least 
one localized state at the phase boundary. For this particular example, numerical diagonalization shows that there are three 
localized states at this boundary. In (e), we take 82+ — IItt/S (orange dot in panel b), so that the quantum walk in all regions 
is characterized by winding number 0. In this case, the probability to find the walker near x = after many steps decays to 0, 
indicating the absence of a localized state at the boundary. 



FigHh- The net result of Eq.(ll) is to make the walker 
hop between sites of a superlattice defined by twice the 
primitive unit cell. The effective Hamiltonian for this 2D 
DTQW takes the form of Eq. (O with the integration 
over k = {kx,ky) taken over the 2D Brillouin zone (BZ) 
of the superlattice. 

We now study the topological properties of the 2D 
DTQWs defined by Eq. pTj) . The corresponding effec- 
tive Hamiltonians lack time-reversal symmetry, and are 
thus contained in the symmetry classes in the bottom row 
of Fig[T] Because U2D is real, this system possesses PHS 
with V = K (see above). With a slight modification, this 
symmetry can be broken and phases in the class labeled 
IQH in FigH] can also be realized (see Supplementary 
Material). These phases are analogous to those of the 
Haldane model j2^, which exhibits an integer quantum 
Hall effect in the absence of a net magnetic field. 

The phases realized by the 2D DTQW, Eq.([ni), are 
characterized by an integer-valued topological invariant 
called the first Chern number. This quantity is defined 
in terms of the unit vector n(k), see Eq.(l5]), as C = 
Tn Jbz'^^^\^ ■ (^fe^'^ ^ ^fc^n)]. Geometrically, the Chern 
number is equal to the number of times n(k) covers the 
unit sphere as k is taken over the 2D Brillouin zone. 
We have numerically calculated Chern numbers for 2D 
DTQWs throughout the full range of spin-rotation angles 



Ol and O2 (see Methods). As shown in Fig|4jD, phases with 
C — Q and with C = ±1 can be realized. 

Similar to the ID case, non-trivial topology in 2D 
DTQWs is manifested in the presence of protected mid- 
gap modes bound to the interface between two topologi- 
cally distinct phases. These gapless modes are analogous 
to the chiral edge modes of quantum Hall systems, and 
are robust against perturbations. To confirm the exis- 
tence of such edge modes, we have used numerical di- 
agonalization to study a non-uniform 2D DTQW on a 
100x100 site triangular lattice with periodic boundary 
conditions, see Fig|4};. We take the spin-rotation angles 
6*1 and 6*2 in Eq. (fTT|) to be site-dependent, with 6*1(2;) = 
6*2(2/) = 37r/2 chosen to realize the C — —1 phase inside 
the red strip 25 < y < 75, and 6*1 (y) = 92{y) = 77r/6 
chosen to realize the trivial C = phase outside. The 
quasi-energy spectrum is plotted in FigHJi as a function 
of the conserved momentum component parallel to the 
interface. Two counter-propagating chiral edge modes 
exist inside the bulk gap. These modes are separately 
localized at the two boundaries between the C = and 
C = — 1 phases, as indicated in FigHj;. 

As described above for ID, these chiral edge modes can 
be probed by performing the spatially inhomogeneous 2D 
DTQW described above with the walker initialized at the 
boundary between two topologically distinct phases. Be- 
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FIG. 4: (a) Translation vectors for the triangular lattice 2D quantum walk defined in EQ. pi|) . (b) Chern number associated 
with the 2D DTQW as a function of the rotation angles 6i and 82. (c) Geometry of an inhomogeneous 2D quantum walk with 
periodic boundary conditions. In the red region, we take 61 = 62 = 37r/2, corresponding to Chern number —1, while in the 
white region we take 61 = 62 = 77r/6, corresponding to Chern number [see colored dots in panel (b)]. Arrows indicate the 
propagation directions of chiral edge modes localized at the two boundaries, (d) Quasi-energy spectrum of the inhomogeneous 
2D quantum walk depicted in panel (c) for a 100 x 100 site lattice. The Brillouin zone for momentum parallel to the interface 
is defined for the doubled unit cell accessed by U2U, Ea. (|lip . Two branches of chiral edge modes connect the upper and lower 
bands. The dotted (solid) line corresponds to the mode localized at the upper (lower) boundary in panel (c). 



cause a general state localized near the phase boundary 
has a non-zero overlap with the chiral edge mode, part of 
the walker's wavepacket will propagate unidirectionally 
along the boundary. Such unidirectional propagation is 
protected by topology, and hence is robust even in the 
presence of an irregularly-shaped boundary. 

Finally, we present a time-reversal-invariant 2D 
DTQW with = — 1, which can realize the quantum 
spin Hall (QSH) phase (see FiglTJ. The realization of this 
phase requires the presence of at least four bands, which 
contain two pairs of time-reversed partners. Therefore, 
we now consider a DTQW where the walker possesses 
four internal states (e.g. a four-level atom, see also exper- 
iment [12] )• We label these four states by a "spin" index 
cr, which takes the values f and J,, and a "flavor" index 
T which takes the values A and B. The time-reversal- 
invariant unitary step operator C/tri is constructed in a 
block-diagonal form, 

^^-^'{"^t Is) W 

where Ua (Ub) only acts on the walker if its flavor index 
is equal to A (B). By fixing Ub — C/J, we ensure that 
Utri is invariant under the TRS operation T = iTyK, 
where Ty is a Pauli matrix which acts on the flavor index. 
As an example, if Ua is chosen according to Eq. pT|) . 
then Ub = i?(-6'i)Tf i?(-02)Tf i?(-6'i)7f ■ Note that 
Tj-^ translates the walker in the direction — (-l-)vi if its 
spin is up (down); i.e. T^^ acts opposite to T^. 

Time-reversal invariant systems in 2D with = — 1 
are characterized by a Z2 topological invariant (middle 
row of right panel in Fig[T]). If 9i and 62 are chosen such 
that Ua is characterized by an odd Chern number, then 
t^TRi realizes a QSH phase with the Z2 invariant equal 



to 1 [19|. Strictly speaking, the effective Hamiltonian 
corresponding to Utri conserves the flavor index r and 
as a result supports topological phases classified by an 
integer Z , rather than the binary invariant Z2. However, 
this additional symmetry can be broken by introducing a 
coupling between A and B states which preserves TRS. 
In this way, the generic Z2 classification can be retrieved 
(see Supplementary Materials). 

DISCUSSION AND SUMMARY 

Because the edge modes bound to interfaces between 
topologically distinct phases in ID and 2D are topologi- 
cally protected, their existence is expected to be robust 
against a broad range of perturbations which may arise 
in real experiments. In particular, their existence is in- 
sensitive to the details of the boundaries, which may be 
sharp or smooth, straight or curved (in 2D), etc. In some 
cases, the topological protection arises from certain sym- 
metries (e.g. chiral symmetry in the ID examples above). 
However, even if these symmetries are slightly broken by 
small errors in the spin-rotation axes and/or angles, the 
edge states are expected to persist due to the absence of 
nearby states inside the bulk energy gap. 

Throughout this work, we have focused on signatures 
of topological phases in single-particle dynamics. How- 
ever, some dramatic manifestations of topological order, 
e.g. charge fractionalization and the quantization of the 
Hall conductivity, appear for specific many-body states 
such as the "filled-band" ground states of fermionic sys- 
tems. To observe these phenomena in DTQWs with mul- 
tiple walkers, analogous many-body states can be pre- 
pared schematically as follows. For special choices of the 
DTQW parameters, the Bloch eigenstates are simple, i.e. 
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local in space, and uniform in spin. By preparing a sin- 
gle filled band comprised of such states, more compli- 
cated filled band states can be obtained through a quasi- 
adiabatic evolution in which the DTQW parameters are 
changed slightly from step to step. Even if an energy gap 
closes along the way, the number of excitations created in 
the process can be controlled by the effective sweep rate. 
In this way, many-body aspects of topological phases may 
also be studied using DTQWs. 

In this paper, we have shown that discrete time quan- 
tum walks provide a unique setting in which to realize 
topological phases in ID and 2D. With only slight mod- 
ifications to the quantum walk protocol which was real- 
ized in recent experiments, the entire "periodic table" of 
topological insulators 2^ 27 [ in one and two dimensions 
can be explored. In addition, we have provided a method 
to detect the presence of topological phases through the 
appearance of robust edge states at boundaries between 
topologically distinct phases. 

Recently, several promising system-specific methods 
have been proposed to realize topological phases using 
cold- atoms |30l- 42 1 . polar molecules [43|, or photons [ii"]. 
Our work advances this emerging field by providing a 
general framework for studying topological phases in a 
wide variety of available experimental systems including 
cold- atoms, trapped ions and photons. By extending this 
work to three dimensions, it may be possible to realize 
new topological phases, such as the Hopf insulator [45], 
which have not yet been explored in condensed matter 
systems. In addition, multi-particle generalizations of 
discrete time quantum walks will open new avenues in 
which to explore the quantum many-body dynamics of 
interacting fermionic or bosonic systems. 



METHODS 

Determination of tlie phase diagram for ID 
split-step DTQW 

The unitary evolution of the ID split-step DTQW, 
Ea.(ITU|. is generated by a Hamiltonian of the form 
of Eq.® with cosE{k) = cos(6'2/2) cos(6'i/2) cos fc - 
sin(6'i/2)sin(6'2/2), and 



nx{k) 
ny{k) 



cos(6'2/2)sin(6'i/2)sinfc 
sin£;(fc) 

sin(6'2/2) cos(6'i/2) + cos(6'2/2) sin(6'i/2) cos k 
sin E{k) 

- cos(e'2/2) cos(6'i/2) sinfc 
sin£;(fc) 



It is straightforward to check that A(0i) = 
(cos(0i/2), 0, sin(6'i/2)) is perpendicular to n(fc) 
for all k. Therefore, the system possesses chiral sym- 
metry ® with r(6'i) = e-*''^(^i) '^/2. As a result, the 



split-step DTQW can be characterized by the winding 
number of n(/c) around the origin, denoted by Z. Using 
the explicit expression for n(fc) in Eq. p^ . we find 
Z = 1 if |tan(6'2/2)/tan(0i/2)| < 1, and Z = if 
I tan(6'2/2)/ tan(6'i/2)| > 1. The spectrum is gapless 
along the lines | tan(6l2/2)/ tan(6'i/2)| = 1. Thus we 
obtain the phase diagram displayed in Fig. [3] (b) . 



Localized states at a phase boundary of 
inhomogeneous split-step ID DTQW 

In addition to the dynamical simulations presented in 
the main text, we have confirmed the existence of topo- 
logically protected edge states with energy E — or 
E — TT in the ID split-step DTQW through an analytical 
calculation for an infinite system with a sharp boundary, 
using 6*2(2;) = 02- for a; < and 92{x) = 62+ for a: > 0. 
Furthermore, we have used numerical diagonalization to 
study the spectrum of a finite (periodic) system on a 
ring which hosts two phase boundaries. In all cases, we 
find that if the phases on the two sides of a boundary 
are topologically distinct, i.e. characterized by different 
winding numbers Z, then a single localized state with 
energy E = or E = n exists at the boundary. 

For smooth boundaries as described in the main text, 
other localized states that are not protected by topol- 
ogy could appear. These bound states always appear in 
pairs with energies E and —E due to chiral symmetry. 
Therefore, when the phases on the two sides of a bound- 
ary are topologically distinct, an odd number of bound 
states appears at the phase boundary ^16||. 



Phase diagram of the 2D DTQW 

Here we briefly describe a general procedure for de- 
termining the phase diagrams of 2D DTQWs. Because 
the value of a quantized topological invariant can only 
change across a phase boundary where a gap closes, we 
first identify the lines in parameter space along which a 
gap vanishes in the quasi-energy spectrum. Once these 
phase boundaries are determined, the topological phases 
between boundaries can be identified by computing the 
topological invariant at any single point within each re- 
gion. For the 2D DTQW with evolution operator given 
by Ea. ([TT|) . we have obtained phase boundaries analyti- 
cally from the spectrum 

cos£'(k) = {cos(02/2) COS01 cos(k • (vi + V2)) 

— sin(02/2) sin 6*1 cos(k • (vi — V2)} cos(v3 • k) 

- cos(6'2/2) sin(k • (vi -I- V2)) sin(v3 • k), 

which gives the lines shown in Fig. |4] (b). We then nu- 
merically evaluated the Chern number within each region 
using C = ^ /g^ d'^k[n ■ {dk^n x dk^n)] with the appro- 
priate expression for n(k). 
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Online Supplementary 
Material 

EXPLICIT DTQW PROTOCOLS FOR ALL 
TOPOLOGICAL CLASSES 

In this section, we provide explicit DTQW protocols 
which can be used to realize topological phases in each 
of the symmetry classes listed in Fig.l of the main text. 
These protocols are summarized in Supplementary FiglS] 
Each DTQW presented in FiglS] can realize both trivial 
and non-trivial phases within a given symmetry class. 
The specific phase which is realized is determined by 
the spin rotation angles which parametrize the quantum 
walk; the system can be driven through a topological 
phase transition by tuning these spin rotation angles. In 
the following, we denote the presence of time-reversal 
symmetry (TRS) with 7^ = ±1 by TRS= ±1, and the 
absence of time-reversal symmetry by TRS= 0. Simi- 
larly, we denote the presence of particle-hole symmetry 
(PHS) with -P^ = ±1 by PHS= ±1, and its absence by 
PIIS= 0. We denote the presence of chiral symmetry 
under the unitary operator F by CII=1, and its absence 
by CII=0. Note that because the chiral symmetry op- 
erator F is unitary, rather than antiunitary, the phase 
of its square does not carry any additional information. 
In particular the transformation F — > e**F results in 
p2 _^ e^'^F^. For ah of the DTQWs considered below, 
the presence of any two of the symmetries {TRS, PHS, 
CH} automatically ensures the presence of the third. For 
example, if a system possesses PHS and CH under the 
operators V and F, then it also possess TRS under the 
operator T = W. 

Doubling Procedure 

Quantum walks with TRS= —1 can be readily con- 
structed from DTQWs with TRS=0 through the "dou- 
bling procedure" used to construct J7TRi[Eq.(12)] in the 
main text. First the walker is endowed with an additional 
two-fold "fiavor" index r which can take either the value 
A or B. We then choose an evolution operator which is 
diagonal in the fiavor index and which satisfies Ub = C/J, 
where C/a(b) is the evolution operator which acts on the 
walker with flavor A{B). With this possibility in mind, 
below we focus on examples with TRS=0. 

One dimensional topological phases 

Symmetry Classes: 

TRS^ 0, PHS= 0, CH= 1 (Z) 
TRS= -1, PHS= -1, CH= 1 (Z) 

The split-step DTQW described by Eq.(lO) of the main 
text realizes the symmetry class with TRS=1, PHS=1, 



CH=1. By an appropriate change of the direction of the 
spin-rotation axes, TRS and PHS can be broken while 
CH is retained. Thus in order to realize the related sym- 
metry class TRS=0, PHS=0, CH=1, we will break the 
PHS of the split-step DTQW. 

In the main text we showed that any DTQW whose 
unitary evolution operator is real possesses PHS with 
V ^ K, where K is the complex conjugation opera- 
tor. The existence of PHS is in fact more general: if the 
two spin-rotations in a split-step DTQW are performed 
around the same axis, and if that axis lies in the xy- 
plane, then the DTQW will have PHS=1. To see this, 
suppose that both rotations are performed around the 
axis (sin(/3, cos(/3, 0). It is then straightforward to check 
that the resulting effective Hamiltonian possesses PHS 
under the operator V = e'^^'^'V/'^Ke'-'^-'f/'^. 

On the other hand, PHS is absent if we choose a ro- 
tation axis that contains a non-zero z component. An 
example of a DTQW with PHS=0 is provided by the 
evolution operator 

KsiOue^) ^T^Ro.{e2)T^Ro.{0i), (13) 

where Ra{0) is a spin-rotation around the axis a — 
-^(0, 1, 1) through the angle 9. Although PHS is absent, 
this system possesses chiral symmetry under the symme- 
try operator F„(6'i) = ie-'^^'-i^^)-''/'^ , where Aa(6'i) = 
(cos(6li/2),^sin(6li/2),^sin(6li/2)). The absence of 
TRS can be checked in the following way. If the energy 
eigenvalues of the two states with momentum k are given 
by ±|£;(fc)|, then TRS= ±1 requires \E{k)\ = \E{-k)\. 
We have explicitly checked that this relation is not satis- 
fied for the DTQW defined by Eq.(IT31), and thus conclude 
that TRS is absent. 

The split-step DTQW above, Eg. p^ . can realize dis- 
tinct topological phases by tuning the spin-rotation angle 
02- For example, the trivial phase with winding number 
Z = is realized with 6i = tt/2 and 62 = 37r/4 and 
the phase with winding number Z = 1 is realized with 
6*1 = 7r/2 and 02 = 7r/4. 

The recent experimental implementation of a DTQW 
with photons [13| employed the rotation operator given 
by the Hadamard gate R-^ — ie^™"-'"'/'^ , with n = 
l/-\/2(l,0, 1). Since the rotation axis contains a non-zero 
z component, we conclude that this "Hadamard walk" 
belongs to the symmetry class TRS=0, PHS=0, CH=1. 

Using the doubling procedure described above, a time- 
reversal symmetric DTQW with TRS=-1, PHS=-1, 
and CH=1 can be constructed based on the DTQW 
defined in Eq. ()13|) . The corresponding evo- 
lution for one step of the DTQW is given by 
d:i&g[Ug{ei,e2), (u£{ei,e2)f]. it is straightforward to 
check that this quantum walk possesses chiral symmetry 
under the operator F = diag[FQ,(0i), F* . By con- 
struction, this DTQW possesses TRS= —1 with T = 
iTyK. Using these two symmetries, we construct a PHS 
operator V = TT with V"^ = -1. 
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FIG. 5: DTQW protocols for each symmetry class of topological phases in ID and 2D. By tuning the rotation angles, all of 
these examples can realize both trivial and non-trivial topological phases within each class. Here T translates the walker to 
the right (left) if its spin is up (down), while TiiTi) translates only the spin up (down) component to the right (left). In 
2D, the translation Ti shifts the walker in the Vi Vi) direction if its spin is up (down), see Fig. 4a of the main text. A spin 
rotation operator Ru{0) rotates the walker's spin through an angle 9 about the axis u £ {y,a,l3}, where a = -^(0, 1, 1), and 
P = (sin(7r/8), cos(7r/8), 0). In most cases, quantum walks with TRS= — 1 and PHS=-1 are obtained by the doubling procedure 
starting from a quantum walk with evolution operator Ua which has TRS=0 and PHS=0. Such cases are separated by dotted 
lines. See main text for descriptions of the relevant symmetry operators. 



Symmetry Classes: 

^5=0, PHS=1, CH=0 (Z-i) 
TRS=-1, PHS=1, CH=1(Z2) 

The construction of a DTQW with TRS=0, PHS=1, 
CH=0 starts from the split-step DTQW with TRS=1, 
PHS=1, CH=1 [see main text, Eq. (10)]. The chiral 
symmetry can be broken by adding extra operations to 
the split-step DTQW. On the other hand, in the previ- 
ous section we showed that PHS can be retained quite 
generally as long as the two rotation axes are the same 
and taken to lie on a;y-plane. 

In order to construct a DTQW with CH=0, we be- 
gin with Eq.(lO) and add an additional spin-dependent 
translation T which translates the walker to the 
right(left) by one lattice site if its spin is up(down) [see 
Eq.(2)]. Explicitly, the evolution operator for one step 
of a representative DTQW from this symmetry class is 
given by 



One dimensional systems with particle-hole symmetry 
exhibit two distinct topological phases 2J| . These two 
phases are indexed by the Berry phase, which can only 
take the quantized values and tt due to the presence of 
PHS. Explicitly, the invariant is given by 

^[-^){^^^,{k)\^k\^l,{k)). (15) 

Here, \i^ib{k)) is the eigenstate in "lower band" with mo- 
mentum k. The DTQW described above can realize both 
topological phases, with the trivial phase {B = 0) real- 
ized for 01 = 7r/2,02 = ■"'/S, and the non-trivial phase 
with B = 1/2 realized for Oi = 7r/2, 02 = 27r/3. 

Using the doubling procedure, we can construct a time- 
reversal invariant DTQW with TRS=-1, PHS=1, CH=1 
based on Eq. pi)) . 



Uss' {91,02) = T^Ryi02)T^Ry{0l)T, 



(14) 



where Ry(0) is a spin-rotation around the y axis through 
an angle [Eq.(3)]. Since Uss' is real, this DTQW retains 
PHS=1 with V = K. The absence of chiral symmetry 
for this walk can be verified by observing that the quan- 
tization axis n(fc) does not lie on a plane which includes 
the origin. Therefore, no single operator F can be found 
which satisfies TH{k) = —H{k)T for all k. 
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Two dimensional topological phases 

Symmetry Classes: 

TRS=0, PHS=1, CH=0 (Z) 
TRS=-1, PHS=1, CH=1(Z2) 

The triangular lattice 2D DTQW defined by Eq.(ll) 
of the main text involves only spin rotations around the 
y-axis. Consequently, the evolution operator [/2D is real 
and possesses PHS=1 with V = K. Therefore, the time- 
reversal invariant DTQW Utri constructed from U2D, 
Eq.(12), is contained in the symmetry class TRS=— 1, 
PHS=1, CH=1. 

As noted in the main text, C/tri is diagonal in the fla- 
vor index t ^ A, B and thus possess an extra symmetry 
related to the conservation of t^. Here we describe a 
more general 2D DTQW with TRS=— 1 which does not 
possess this additional symmetry. The operator for one 
step of this modified time-reversal invariant DTQW is 
given by 

where Ua{Ub) acts on the walker if its flavor index is 
A{B). The rotation e"*'^"'^'' '^/^ exphcitly introduces mix- 
ing between the A and B flavors, and thus breaks the 
conservation of r^. 

This DTQW is characterized by TRS=-1 with the 
symmetry operator T = iTyK if Ub is chosen according 
to Ub = C^J- If the Chern number associated with Ua 
is odd, then J7tri' with ip = realizes a non-trivial QSH 
topological phase. Because this phase is protected by 
TRS, the presence of a small > can not take the 
system out of this phase. 

Symmetry Classes: 

^5=0, PHS^Q, (Z) 
TRS^-1, PHS^O, CH^O (Z2) 

The existence of topological phases characterized by a 
non-zero Chern number does not rely on the presence of 
PHS. Therefore, the topological phase with Chern num- 
ber 1 in the TRS=0, PHS=1, CS=0 symmetry class can 
be directly transformed to the corresponding phase in 
the TRS=0, PHS=0, CS=0 symmetry class by a pertur- 
bation which breaks PHS. Such a perturbation can be 
achieved by changing the rotation axis for the second ro- 
tation stage in Eq.(ll). The resulting DTQW single-step 
evolution operator is given by 

U^c{di,02) = TsR{ei)T2Rp{02)TiR{ei), (17) 

where i? is a spin rotation around the y axis, and Rp{0) 
is a spin rotation around the axis (3 = (sint/s, cos(^, 0) 
with <p = 7r/16. The operators {T^} correspond to 



spin-dependent translations along the directions {v^}, 
as defined in Fig. 4a. The absence of PHS is confirmed 
by examining the relationship between energy eigenval- 
ues |£'(k)| and |i?(— k)|. The presence of PHS implies 
|£;(k)| = |.E(-k)|. This condition is violated for DTQW 
([T7|. Therefore this system does not possess PHS. 

DTQW P7)) realizes both topologically trivial and non- 
trivial phases with zero and nonzero Chern numbers. 
For example, the choice Oi = 62 = 37r/2 generates the 
phase with Chern number —1, while 9i = 62 = 77r/6 
corresponds to the phase with Chern number 0. Since 
PHS is absent, this DTQW belongs to the class with 
TRS=0, PHS=0, CS=0. The related time-reversal in- 
variant DTQW constructed by applying the doubling 
procedure to this walk has TRS= -1, PHS=0, CH=0. 



Symmetry Class: 

TRS=0, PHS=-1, CS^O (Z) 

Quantum walks with PHS=— 1 can be constructed 
through a doubling procedure similar to that used to 
construct DTQWs with TRS=-1. Consider the block- 
diagonal evolution operator 

c/p,, -{\^ I). m 

where Ua{Ub) acts on the walker if its flavor index is 
A{B). If we choose Ub^U\, then the resuhing DTQWs 
possess PHS=— 1 with V — iTyK. By choosing Ua ac- 
cording to Eq. ()17|) with parameters to give a Chern num- 
ber of 1, Eq. ([T5| produces a DTQW which realizes 
a non-trivial topological phase in the symmetry class 
TRS=0, PHS=-1, CH=0. 



REALIZATION OF 2D TOPOLOGICAL PHASES 
ON A SQUARE LATTICE 

In the main text, we have provided examples of 2D 
DTQWs which realize topological phases on a triangular 
lattice. However, these DTQWs can also be implemented 
on a square lattice, as we explain below. A square lattice 
may be easier to realize in some experimental implemen- 
tations, such as cold atoms in optical lattices. 

Generally speaking, the phase diagram of a DTQW is 
determined by the amplitude for the walker to hop from 
one site to another after one complete step of the evolu- 
tion. Thus, as long as the hopping amplitudes between all 
pairs of sites are preserved, geometrical deformations of 
the lattice do not change the phase diagram. Therefore, 
the phase diagram of a DTQW is insensitive to geometric 
deformations of its host lattice. 

In particular, 2D DTQWs with non-zero Chern num- 
ber can be realized on a square lattice by replacing the 
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— ► spin up 
— ► spin down 



+W2a 
-W3 




-W2 





FIG. 6: Translation vectors for 2D DTQW on a square lattice 
that realizes the phase diagram of Fig. 4 in the main text. 
Crucially, V3 satisfies the relation V3 = vi — V2 just as in the 
triangular lattice realization. 



translations along the vectors {v^} on the triangular lat- 
tice in Eq.(ll) with the vectors {w^} shown in Fig|6l 
Here, wi = (1, 1), wa = (0, 1) and W3 = (1,0). This pro- 
tocol is obtained simply by "shearing" the lattice used in 
Eq.(ll) and Fig. 4a. Note that the diagonal translation 
along Wi can be implemented by a compound translation 
along (1,0) followed by a translation along (0, 1). 



